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I. INTRODUCTION 



In a previous work we have treated the problem of a rigorous determination of the upper 
bound for the kinetic density functional of a system of electrons [lj] . There, starting from 
the most general expression of a many-body wavefunction in polar form, a final formula 
was obtained where the kinetic functional consisted of two terms. The first, linked to the 
phase factor of the wavefunction, we have shown, can be bounded from above by a term 
proportional to the Thomas-Fermi one, the second corresponds to the Weizsacker term. A 
third term, corresponding to what is known as the non-local Fisher information functional, 
was neglected since, being quantitatively much smaller than the other terms it was not 
relevant for the discussion carried on there. However, although quantitatively not relevant, 
such a term plays a fundamental role regarding important structural properties, such as 
the atomic shell structure, and it is strictly related to the electronic correlations implicit in 
the kinetic functional [2I (for this reason we have renamed it as "kinetic correlation term"). 
The difficulty in treating it is due to the fact that an explicit rigorous expression does 
not exists and to find a suitable one represents a challenging task [4(. To our knowledge, 
only few attempts have been made in this direction and they were restricted to highly 
simplified situations |5j. Clearly, for a detailed electronic description of condensed matter 
systems, this terms may play a key role. For this reason we focus the attention on such 
a problem and propose a general procedure, physically well founded and mathematically 
rigorous, to obtain an analytic expression. The practical importance of this problems stems 
from the fact that, as it is discussed in 1], there is a renewed interest towards energy 
functionals which are nowadays used in flexible quantum (based) computational methods, 
such as the linear-scaling real-space kinetic energy functional method, where the kinetic 
energy is calculated as a functional of the electron density. Such a method is also known 
as Orbital Free Density Functional (OFDFT) (see e.g. 6|). These modern methods allow 
in principle for the treatment of much larger systems and longer time scales than those 
based on standard Density Functional Theory in the Kohn-Sham formulation 0]. The 
crucial point of OFDFT, is that it is in principle exact, provided that the kinetic functional 
is rigorous. This means that the designing of rigorous kinetic density functionals would 
allow to explore the full power of the Hohenberg-Kohn theorem for extended systems. In 
the light of the discussion above, the search for a rigorous kinetic functional is not simply 
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interesting for its own theoretical implications and understanding of fundamental aspect of 
quantum mechanics, but also for its large importance in practical calculations in condensed 
matter. The paper is organized as follows: the many-body wavefunction approach to obtain 
a 3iV-dimensional functional (with N, the number of electrons), which was presented in 
our previous work, is reviewed and extended to the case where the electrons are explicitly 
interacting. Next the reduction to a 3-dimensional case for each term appearing in the 3N- 
dimensional expression is treated. While the term related to the wavefunction phase factor 
and the Coulomb interaction are known, the part related to the kinetic correlation is new 
and represent the core of this work. 

II. THE iV-PARTICLE WAVEFUNCTION APPROACH FOR THE HAMILTO- 
NIAN PROBLEM 

Let us consider the many body Hamiltonian operator for a system of N interacting 
electrons: 

1 N 1 1 

t=l j>i I 1 J I 

where for simplicity the Planck constant h, the electron mass m and the electron charge e 
are taken to be 1. Let us consider the most general form of the wavefunction solution of the 
stationary Schrodinger equation for an arbitrary state of energy e: 

Hip = eip (2) 

Such a general form can be written as: 

V(ri, ,r N ) = e iS <* ^9( ri , , r^) (3) 

with S and O real functions in an ^-dimensional space of volume Q = Ylf =l Ui = u N , 
with uj being the single particle volume. We also require S to be symmetric and to be 
antisymmetric with respect to any particle pair permutation, tp satisfies the normalization 
condition: J n ifj*ifjdr 1 ...,drN = 1. The spin variables are not explicitly considered, however, 
their effects will be properly taken into account and specifically mentioned when it will be 
required. Our aim is to use Eqj2] and the form and properties of ip to derive an energy 
functional E[p] = J F[p(r)]dr where _F[p(r)] is the energy density written as a functional of 
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the one particle electron density, p(r), denned as: 

p(t)=p(t!) = N tp*(r 1 ,....,r N )tfj(r 1 ,....,r N )dr 2 dr N (4) 

J ffiJV— 1 

where Qn-i is the A" — 1-dimensional volume of A" — 1 particles, and / p(r)dr = N; since 
the particles are indistinguishable, the index 1 can be exchanged with any other of the 
remaining N — 1 particle index, for the rest of the work this formalism will be kept. E[p] is 
the equivalent of e in Eqj2] written in terms of p(r); thus we should obtain: 

= J {$*Hil)) dri.... dr N = e = E[p] = J F[p(r)}dr. (5) 



E\ 



h = \H e 2 (r u ..,v N )\V t S(r u ..,v N )\ 2 + l £ | V^, .., r N ) | 2 + 

1 i=l,N 1 i=l,N 

i N l,N ft 2( r „ \ i 

+ .' 'I": i -I E V?8'(r 1 ,..,r J ,). (6) 

Eqj6] is defined in a 3iV-dimensional space; in order to derive F[p(r)}, we shall reduce the 
problem to a 3-dimensional space. Note that the last term on the r.h.s. vanishes after 
integration, thus it will not be relevant in this context. The first step to reduce the problem 

to a 3-dimensional space, is to redefine 2 (r 1; r^v) in terms of a one particle probability 

density and an (N — l)-particle conditional probability density: 

e 2 (r 1 ,....,r^) = ^/(r 2 , (7) 

where = ^j^- is the one particle probability density and /(r 2 , , rjv/rj) is the N — 1 

electron conditional (w.r.t. r x ) probability density, i.e. the probability density of finding an 

N — 1 electron configuration, C(r 2 , , r N ), for a given fixed value of rj. The function / 

satisfies the following property: 

/(r 2 , , r N /T 1 )dr 2 dr N = l\fr 1 (8) 

iV-l 

as before, due to the electron indistinguishability, the index 1 was chosen arbitrarily. Another 
equivalent factorization of G 2 , which it will be used for simplifying the electron-electron 
Coulomb term, consists of writing 2 as a product of a two particle probability density: 



POi,r 2 ) 



Af Jn N _2 



/ ^*(r 1 ,r 2 ,r 3 ....,r Ar )V'(ri,r 2 ,r3....,r 7V )rfr3 dr N (9) 



equivalent to the diagonal element of the second-order spinless density matrix 7], and an 
N — 2 particle conditional probability density: 

/(r 3 , ,Wn,r 2 ). (10) 

Equivalently to the factorization shown before, in this case we have (see e.g. jsj): 

J p(r 1 ,r 2 )dr 1 dr 2 = N (11) 



and: 



/ /(r 3 , , rjv/ri, r 2 )dr 3 ....dr N = lVri, r 2 . 



(12) 



As a consequence of the definitions given above, the following identity is satisfied: 



N 



—f(r 2 ,..,r N /rx) 



N 



/(r 3 ,..,r Ar /r 1 ,r 2 ). 



(13) 



Once again, since the particles are indistinguishable, the indices 1 and 2 for the one and 
two particle density were chosen arbitrarily and any other choice would be equivalent. In- 
serting EqjT] and the r.h.s. of Eq{13] into EqEJ using the properties of / shown above, and 
remembering that the particles are indistinguishable, we obtain (for the kinetic terms see 
also HQ): 



/ (ij>*Hil>)dr 1 ...dT N = ± [ J2® 2 (^-^N)\V l S(T 1 ,....,r N )\ 2 dv 1 ...dr N + ]- f |V/ ' (ri 



1 P(ri) 


/ 




Jn 



\V ri f(r 2 ,....,r N /r 1 ) 

/(r 2 , ,W r i) 

, (N-l) 



P( r i) 
-dr 2 ....dr 

P(*"l,r 2 ) 



dridr^U) 



2 | i*i — r 2 

The last term on the r.h.s. of EqJTH is obtained in the following way: first we use the 
property of particle indistinguishability and rewrite the Coulomb term in Eqj6] as : 



if ff mn, 



-dri...G?rjv 



iV(iV-l) /■ |e( ri ,....,r^ 



ri - r 2 



-dri...dr 



N 



(15) 



where the indices 1 and 2 are, as usual, chosen arbitrarily, then we use the decomposition 
of 6 2 as a product of the two particle density p(ri,r 2 ), and the two particle conditional 
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probability density, /(r 3 , rjv/ri, r 2 ), and obtain: 

N(N-l) , |e (ri r w )|» (W-l) , /^ dridr2 



n |ri — r 2 | z Ju; |ri — r 2 

/ /(r 3 ,....,rjy/ri,r 2 )tir3 rfrjy (16) 

EqlTBI can be simplified to a two particle expression, since 

/ /(r 3 , ,r N /ri,T 2 )dr 3 dr N = lVr l5 r 2 . (17) 

so that we have: 

JVfAT-1) , |e(r, r„)|« = (JV-1) , , r<JulA dridr2 (18) 

2 |ri — r 2 | 2 JwJw|ri — r 2 | 

III. FROM THE iV-PARTICLE TO THE ONE PARTICLE REPRESENTATION 

As one can see from EqJT^l there are two kinetic terms which are written in an iV-particle 
representation, namely, the phase factor term: 

1 N 

/ -Y l ^(T 1 ,...,T N )\V i S{T 1 ,....,T N )\ 2 (19) 

and the correlation term also known as the Fisher non-local information functional 



ii 
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1 f ( x , f |V ri /(r 2 , ^iv/rx)! 2 

- / p(ri)dri / r dr 2 dr N (20) 

o Jui Jn N _ 1 j [T2, ,rjvj 

and the aim of this work is to reduce them to a one-particle representation; the first term 

was already considered in a previous work [1] and the results will be reported here for 

completeness; the second term deserves a particular attention, since to our knowledge only 

very few attempts, restricted to specific and simple cases 5|, were made to give an explicit 

formula; we will propose a rather general procedure to obtain an analytic expression. The 

third term, the electron-electron Coulomb interaction, which, as we have seen before, is a 

two particle term, and, generalizing the index 1 and 2, can be written as: 

E „ = ^H 1 1 ^l dld , (21) 

2 Jui Jui' |r — r'| 

it has been extensively treated in literature and the related results will play an important 
role in the procedure of reducing the iV-particle expression of Eqj20] to a single particle 
functional expression. 



A. Reduction of J n £)t=i n 2 |Vi/S , | 2 dri...rjv to a one particle representation 

In our previous work [l| we have shown that in the limit of a slowly varying density we 
can reasonably conjecture the following inequality: 

|ViS| < |PfJ (22) 

which means that the momentum of the i-th particle is less (or equal) to the Fermi mo- 
mentum Pp i assigned to the i-th particle (i.e. the maximum value of the momentum an 
electron, under such conditions, can have). We have shown that Eqj22] leads to the following 
inequality for the kinetic term related to the phase factor S: 

. N 

/ ^e\T 1 ,...,T N )\V i S(r 1 ,...,T N )\ 2 dr 1 ...dr N <C / p§(r)dr. (23) 

J Q ■ ^ J LJ 

In %\2/3 

Where C = - — ^ — . We have also discussed the reasonability of our results beyond the 
slowly varying density regime and how our upper bound compares to others available in 
the literature. The effects of the spin are implicitly incorporated into the definition of P^, 
where the spin degeneracy enters as a scaling factor. On the other hand, any possible spin 
(or spatial) correlation effect implicit in S can be neglected, because the energy assigned to 
each electron in the upper bound term is so large, or equivalently the overestimation of the 
energy is large enough, that any possible correlation energy can be reasonably neglected in 
Eqj23l In this case the term "reduction" means that the 3iV-dimensional term is reduced to a 
3-dimensional term, expressed through the quantity p(r), which bounds the true expression 
from above. This is the first step in reducing E[ip] = f n ip*H^dri...dr N to a reasonable 
density functional. 

B. The electron-electron interaction 

The electron-electron interaction energy E ee = ( - JV ~ 1 - ) / w f^, ^-£'1 drdr', can be separated 
in two parts, the (classical ) Hartree term: 

E„=ff ^Idrdr' (24) 

and the exchange and correlation term: 

E xc = [ p{v)e xc {v)dr (25) 
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where e xc (r) is the exchange and correlation energy density (see e.g. |7|). In order to 
describe the nature of the exchange and correlation term in a simple and compact way, 
we refer to a series of seminal papers by Wigner and his coworker Seitz [13I . 11^ . In their 
work on the interaction of electrons in metals, they treat the case of uniformly distributed 
interacting electrons. They argue that the interaction energy of electrons with parallel spin is 
determined by the classical Coulomb interaction due to the space charge plus the exchange 
term. This latter, also known as Fermi correlation [151], is the consequence of the Pauli 
exclusion principle, and expresses the fact that two electrons with the same spin stay away 
from each other; such an effect is also called a Fermi hole. In contrast, the interaction of 
electrons with antiparallel spin contains the classical Coulomb term, plus the correlation 
term. The latter arises from the fact that electrons with antiparallel spin tend to condense 
in the same place but the electrostatic repulsion keeps them apart. We should also add 
that the exchange and correlation effects contained in the electron-electron term, are, by 
now, well understood. Although their satisfactory description is by far not easy achievable, 
there is a rich literature dealing with reasonable descriptions (numerical, analytical, semi- 
analytical) of such a quantity (see e.g. articles in [lij and references therein). It is not the 
purpose of this work to propose anything new regarding this part; on the contrary, it is 
the availability of various analytical formulas for the exchange and correlation term which 
makes it possible to formally derive results on the kinetic correlation term. 



C. The non-local Fisher Information Functional or Kinetic Correlation term 



After having treated the kinetic term related to the wavefunction phase factor and having 
given a general description of the direct electron-electron interaction, we now come to the 
core of this work; to give a general prescription for finding an analytical functional form for 
the non local Fisher information functional or kinetic correlation term: 

|V r /(r 2 ,.... r7V /r)' 2 



I [ p(r)I(r)dr = - f p(r) 

O Jul O JUJ 



-dr 2 ...dr N 



dr. 



(26) 



iq n _x f(r 2 ,....r N /r) 

The main point is to give physically as well as mathematically well founded conditions which, 
when satisfied, lead to a valid expression for I(r). This term is qualitatively very important 
since it is responsible of the atomic shell structure 0, 5], thus to have a valid description 
of it automatically means to have an accurate description of the electronic structure. Inter- 
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estingly, Ludena and Karasiev[3[ intuitively link this term to the exchange and correlation 
energy, although a prescription for how to find this link is not given. We start from the 
same intuitive idea and develop it further. The essence of the idea developed here is that a 
gas of electrons, as long as these are not related by an explicit interaction, does not display 
correlation effects, except those of the basic Fermi distribution, or occupation of levels, for 
non interacting electrons. However, when a direct electron-electron interaction is switched 
on, then the electrons start to feel the presence of all the others in a non trivial way. This 
means that their spatial distribution is heavily governed by the fact that they try to opti- 
mize such correlations (interactions) in a way that the optimal spatial arrangement fulfills 
the minimum energy requirements. For this reason the key quantity determining the many 
particle electron distribution is the electron-electron interaction: 

^ ee (r,r') = f^. (27) 
|r — r | 

The expression of Eqj27]is a two-particle function, averaged over the rest of iV — 2 particles, 
where the effects of these latter on the electrons at r and r' is implicitly described. As we 
have seen before, although we do not have an exact expression for E ee (r,r > ), nevertheless 
several valid approximations are given in an analytic form. Conversely, very little can 
be found in the literature for I(r) due to the difficulty in writing an explicit expression 
for /(r 2 , rjv/r) [4(. For this reason, since we cannot write directly /(r 2 , rjv/r) in an 
explicit way, we can think of transforming /(r2, rjv/r) into some function of E ee (r, r'); 
this approach, although indirect, can lead to a physically well-founded analytic expression. 
As said before, /(r 2 , rjv/r) is the probability density of finding N — l particles in a certain 
configuration once one particle r is fixed; it is reasonable to think that it must be linked 
to the electron-electron interaction of EqfSTJ For the reasons given above we propose the 
following conjecture: 

/(r 2 ,...,r iV /r) = UZ 2 hi(E ee (r, r*)) (28) 
where hi(E ee (r,Ti)) must satisfies the following conditions: 

(a) f u hi(E ee (r, r^dr* = 1; Vr 

(b) h i (E ee (r,r)) = (29) 

(c) UgMEeeir,^)) < l;Vr,r, . 
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Condition (a) comes directly from EqJHJ condition (b) is the direct consequence of the anti- 
symmetry of ip, while the last condition expresses the fact that /(r2, rjv/r) is a probabil- 
ity density, thus cannot be larger than one. The function H^L 2 hi{E ee (r, r,)) in this context 
should be interpreted as the probability density which expresses how the particles dis- 
tribute themselves w.r.t. the r-particle (fixed), while experiencing each other (except r) in 
an average way implicitly contained in E ee (r,r') through p(r). The conditions in Eqj29l are 
general, and can be used to design possible functional forms for I(r). In the next section we 
propose some functional forms and discuss the consequences. 

1 . The exponential form 

In this section we will treat the case when the function hi = hi(E ee (r, r^)) has the form of 
an exponential. By this ansatz and its related mathematical prescription we will gain some 
interesting insights into, possibly, physical effects related to the kinetic correlation term we 
are treating. The functional form we propose for / has the following expression: 

TlgMEeeir, r,)) = n^-^e-^O^) (30) 

where 



e -E(r) = / e -See(r,r l)rfr ._ (31) 

J Ul 

This function naturally satisfies the requirement (a) and (c) of Eqj29j in fact 

e^ (r) x / e-^'^dr, = l;Vr,i. (32) 

J UJ 

and as a consequence, because of the particles' indistinguishability, we have: 

/ Ut 2 e {N - mr) e- E - ir ^dr 2 ...dv N = 1; Vr. (33) 

Condition (c) is satisfied, since e E ^ is the automatic normalization factor, and thus: 

-_B ee (r,rj) 

=7T— <l;Vr,r< (34) 

e -E(r) — ' ' 1 V ' 

and again, this can be straightforwardly generalized to the complete T\.f =2 hi(E ee (r, r^)) func- 
tion. We will come back to condition (b) of Eql29] later on, for the moment let us jus- 
tify why the exponential form should be physically sound. The starting point is that the 



10 



electron-electron interaction is a positive definite function, E ee (r,r') > 0; Vr, r'. The phys- 
ical principle we use is that the electrons try to distribute themselves in a way that their 
spatial arrangement is consistent with the optimal electron-electron total interaction energy. 
In this context the conditional probability density, /(r 2 , ...rjy/r) = TlfL 2 z^z^w^i expresses 
the fact that for a configuration of high interaction the probability is low because electrons 
tend to avoid each other; on the other hand, for a lower interaction energy the probability 
is higher, because the electrons can distribute in a more optimal way. So far, on the basis 
of the physical intuition and mathematical requirements, our ansatz seems rather reason- 
able; however in this form, /(r2, ...rjv/r) does not satisfy a major symmetry requirement, 
namely condition (b) of Eql29l In fact 0(ri, ..rjv) being antisymmetric under any particle 
pair permutation, it would require that 2 (r, r 2 , ...r, ...rjv) = 0. We should remember that 
2 (r, r 2 , ...r, ...ttv) = p(r)/(r 2 , ...r N /r), which implies that /(r 2 , ., r, .., r N /r) = 0; this is 
impossible in the form we have expressed /(r 2 , ...r^/r). To understand this point we have 
to look at the explicit expression of E ee : 

E ^') = y^7\ ( 35 ) 

thus: 

lim£ee(r,r') = (36) 

r^r' 

i.e. the probability of finding two electrons in the same position is zero, and thus the related 



interaction is zero too. This can be also seen, for example, in Ref.[15[], where E ee (r, r') is 
explicitly written in term of classical Coulomb interaction plus the exchange and correlation 
term; for r — > r' the exchange and correlation part becomes equal to minus the classical 
Coulomb term, and thus they cancel each other. For our expression of /, instead, we 
must require E ee (r,r') — > oo for r' — > r so that /(r 2 , .,r, ..,r N /r) = 0, however this is not 
possible because E ee (r, r) = thus we cannot use E central quantity. There is 

also another necessary condition which must be satisfied by /(r 2 , .., rjv/r), and it involves 
directly the kinetic correlation term; in fact, this latter must be finite whenever any one of 
the Tj electrons approaches the r-electron; i.e. 

|V r /(r 2 ,..,...r, ...,r7v/r)| 2 

77 7^ = fmite (37) 

/(r 2 , ..,r, ...,1-jv/r) 

where finite is in general either zero or a finite positive number or even a non diverging 
function of the other variables than r. Eqj37] implies that the numerator should go to zero 
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with the same order of, or faster than, the denominator. This mathematical requirement 
together with condition (b) of Eqj29] suggest that instead of the full electron-electron inter- 
action, E ee (r, r'), we should use only the classical Coulomb term, also known as the Hartree 
term, £#(r, r') = p j^_ p ^ - . It can be easily verified that it satisfies the conditions (a), (b), 
(c) of Eql29l and leads to a finite value (zero) for the kinetic correlation term in the limit of 
i"j — > r. In fact in this case one has: 

2 



|V r /(r 2 ,.., rjV /r)| 2 
/(r 2 ,..,r7v/r) 

The r.h.s. of Eql38l becomes: 

e (N-l)E(r) u N_^ e -E(r,r z ) 

and thus we have: 

lim e (N-l)E(r) U N -Eir,^) . 



jjJV e (JV-l)E(r) e -B H (r,ri) 



n,C2e (JV_1)B(r) e" %(r,r<) 



(38) 



N 



(N -l)V T E(r) -^VrEH^n 



i=2 



(39) 



N 



(N - l)V r E(r) - "£V r E H {r, ri 



i=2 



0;Viv 



(40) 



The limit in EqH0]is zero because while (N — l)Vr-E'(r) — Y^f=2 V T Eg(r, r$) diverges for 
r, — > r, in the worst | n with n a finite positive integer, the exponential 

e -E H {r,Vi) tends to zero in a much faster way than the divergent term. The fact that I(r), 
written using the exponential form of /, does not diverge for — > r, together with the 
other mathematical properties discussed before, show that indeed the exponential ansatz is 
mathematically and physically sound, thus we can proceed to the core of this work, that is 
to write an explicit formula for I(r). From Eqj38l and Eqj39]we obtain: 

2 



I(r) = (N-iy 



V P Mr 



N 



-2(#-i)E 



i=2 



e (.N-2)E(r) Ui ^ e -E H (r,r t ) dTi 



N 



,E(r)-E H (T,Ti) 



V r E(r) ■ V r E H (r,r 



i=2 



e (iV-2)S(r) nj/ . e -BH(r,r l ) drj 



e £ W-^( r ' r ^|V r ^(r,r i )| 2 rfr i + 



+ EE 



e B(r)-£(r,r,) e £(r)-J?(r,r J ) x 



x \V t E{t, n) ■ V r £(r, Tj )\ dr^41) 



Using the indistinguishability, taking into account that all the expressions of the form 
/ e E ( r ^~ E ( r ' r ^dri = 1; Vr and factorizing all the r^'s that do not contain the variable r$ present 
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in V r -E(r, rj), we obtain: 

J(r) = (N - l) 2 £(r) 2 - 2(iV - l) 2 / e S(r) - £(r ' r,) V r £(r) • V r £(r, r')dr' + 

+ (N-1) ( e^ (r) - £(r ' r ' ) |V r ^(r,r , )| 2 rfr / + 
+{N - 1){N - 2) / / e S(r) - £(r ' r ' ) e S(r) - £(r ' r " ) V rJ B(r, r') • V r £(r, r")dr'dr". (42) 



The last term on the r.h.s. of EqHH characterized by a sort of double independent non 
locality, (r' and r"), can be reduced to a simple non locality since: 

V r £(r, r') ■ V r £(r, r") = d x E(r, r')d x E(r, r") + d y E(r, r')d y E(r, r") + «9 2 £(r, r')d z £(r, r") 

(43) 

thus: 

(N - 1){N - 2) / / e^ r )- £ ( r ' r ' ) e^ r )- £ ( r ' r " ) V rJ E(r, r') • V r £(r, r")dr'dr" = 

= (N — 1)(N — 2)e 2E{r) ]T ! e E{ry) d Xk E(r, r')dr' f e E ^d Xk E(r, v")dr" (44) 

k=l,3 Ju> Ju 

where x± = x; x 2 = y; x 3 = z, and from Eqj4"4"j since d Xk E(r, r') and d Xk E(r, r") are indepen- 
dent from each other, one obtains: 

e 2E(r) J- f e E ^d Xk E(r,r')dr' [ e E ^d Xk E(r, r")dr" = e 2E ^ £ If e E ^d Xk E(v,r')dv' 

k=l,3 Juj Ju) fc=l,3 ^ w 

(45) 

Inserting EqSS] in Eqj4"2"l we obtain an analytic expression for /(r), since all the terms 
involved in EqJ42] and Eqf45] are known explicitly: 

p(r)p(r') 



and; 



E„(t,t') 



E(t) = -log 



r — r 



l e -E H (r,r>) dr , 

J UJ 



(46) 



(47) 



The physical interpretation of the results obtained is the following, the non local Fisher 
information functional or kinetic correlation functional, can be well described by the classical 
Coulomb interaction between the electrons. Such a conclusion is based on the fact that the 
exponential form is physically "reasonable", however the related mathematical prescriptions 
do exclude that the exchange and correlation interaction plays a role in this case. Rather, it 
leads to the conclusion that the classical Coulomb interaction naturally dominates the kinetic 
correlation effects since it perfectly matches the physical arguments and the mathematical 
requirements. 
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D. Discussion and Conclusions 



The proposed exponential form for the conditional probability has been obtained on 
the basis of well founded physical arguments and rigorous mathematical prescriptions. In- 
terestingly, starting from a rather different point of view and aiming at a different goal, 



Ceperly 17] and Jastrow [18j], choose, on the basis of physical intuition, the correlation part 
of a fermionic wavefunction as e~ u ^ ri ~ Tj ^ , i.e. a repulsive pseudopotential. Moreover, Wigner 
and Seitz [14J], and Wigner [ijj], conclude that the conditional probability part contained 
into a fermionic wavefunction, must be of the form /(|r^ — r^|). This suggest that indeed, 
the physical principle we have used, that is a high Coulomb interaction corresponds to a 
less favorable situation compared to the situation of a weaker interaction, is a justified argu- 
ment. Moreover, the exponential form given to the conditional probability, which resembles 
the behavior just explained and fulfills all the mathematical requirements, although not 
exact, cannot be far from the exact one , and thus, at least qualitatively, also the func- 
tional obtained is reasonably rigorous. A practical advantage of this approach is that it is 
not difficult, in principle, to implement this idea in a self-consistent calculation for atoms 
with a statitically relevant number of electrons (e.g. metals), where the Thomas- Fermi and 
Weiszacker terms describes well the the local part of the kinetic energy. In this case one 
can use the spherical symmetry and then compare the results with those obtained using, 
for example, the Hartree-Fock method. This test will certainly show possible advantages or 
disadvantages of this approach and thus inspire new ideas to improve it. 

Acknowledgments 

We thank M.Deserno, R.A.Mosna and I.P.Hamilton for helpful comments. 



[1] L.Delle Site,J.Phys.A 38, 7893 (2005). 

[2] E.V. Ludena, V.V.Karasiev and P.Nieto, Theor.Chem.Acc.110, 395 (2003). 

[3] E.V. Ludena and V.Karasiev, Kinetic Energy Functionals: History, Challenges and Prospects, 
in: K.D. Sen (ed.) Reviews of Modern Quantum Chemistry. A Celebration of the Contri- 
butions of Robert G. Parr, World Scientific, New Jersey, 612, (2002). 



14 



[4] M.Kohout, Int.J.Quant.Chem. 87, 12 (2002). 

[5] J.L.Gasquez and E.V. Luderia, Chem.Phys.Lett. 83, 145 (1981). 

[6] S.Watson and E.Carter, Comp.Phys.Comm. 128, 67 (2000). 

[7] W.Yang and R.G.Parr, Density Functional Theory of Atoms and Molecules, Oxford University 
Press, New York, 1989. 

[8] P.W.Ayers, J.Math.Phys. 46, 062107 (2005). 

[9] R.A.Fisher, Proc. Cambridge Philoso.Soc. 22, 700 (1925). 

[10] S.B.Sears, R.G.Parr and U.Dinur, Isr.J.Chem. 19, 165 (1980). 

[11] A.Nagy, J.Chem.Phys. 119, 9401 (2003). 

[12] E.Romera and J.S.Dehesa, J.Chem.Phys. 120, 8906 (2004). 

[13] E.Wigner, Phys.Rev. 46, 1002 (1934). 

[14] E.Wigner and F.Seitz, Phys.Rev. 46, 509 (1934). 

[15] S.K.Gosh and B.M.Deb, Phys.Rep.92, 1 (1982). 

[16] A Primer in Density Functional Theory, C.Fiolhais, F.Nogueira and M. Marques (Eds.), Lec- 
ture Notes in Physics, 620, Springer- Verlag, Berlin Heidelberg 2003. 

[17] D.Ceperly, Phys.Rev.B 18, 3126 (1978). 

[18] R.Jastrow, Phys.Rev. 98, 1479 (1955). 



15 



